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Local unitary invariants of multipartite states fall into families related by the tracing-out of 
subsystems. In the case of pure qubit systems, there is a family that accounts for about half the 
total number of invariants and is closely connected to multipartite separability. One way to define 
this family is to give pure states the structure of an algebra, and define a log function in this algebra. 
The coefficients of the Taylor expansion of this log function, which are polynomials in the coefficients 
of the states, are cumulants. When twirled by local unitaries, these yield invariants. The traditional 
cumulant, which is a function of random variables, vanishes if its arguments belong to two or more 
independent sets. The equivalent of this in our context is that certain cumulant-invariants vanish 
when a state is separable. 



I. INTRODUCTION 

Given a multipartite quantum state, the action of local unitaries maps out an orbit of locally equivalent states. 
We can think of such orbits as forming an entanglement space, £. The orbits can in principle be characterised by a 
complete set of invariants of this group action, and the problem of finding complete, algebraically independent sets of 
invariants has been much studied [1, Il5l - [l7i [20 . [23j . The picture that emerges is not altogether encouraging because 
the number of invariants grows exponentially with the number of systems. However, there is a redeeming feature: an 
invariant in n systems gives rise to a set of invariants in n + 1 systems via a tracing-out operation. This means that 
any invariant generates a family that grows exponentially with the number of systems, and one can therefore hope to 
grasp the structure of the space through these families. 

I illustrate this by defining a family of invariants for pure qubit states that grows exponentially with n and asymp- 
totically accounts for half the total number of invariants. The key to the construction of this family is the set of joint 
cumulants @, 0, EH HH, H^, [26|. Cumulants are most commonly encountered as statistical tools, or as ingredients 
in cluster expansions. For our purposes, they are simply polynomials in the coefficients of states that have certain 
desirable properties; for instance, they are closely related to the separability of multiparty states. Cumulants can be 
introduced in an attractive if unorthodox way by giving the space of vectors in (C d )® n the structure of an algebra in 
which a Taylor series and hence analystic functions can be defined. In particular, one can define a log function with 
the property that 

log(^>®|^) = log(|V»+log(|0}), 

for all \ip) and \4>); see © and (|TU)) for a more precise statement. The coefficients in the log-expansion are cumulants, 
which is how the connection between cumulants and separability comes about. 

This algebra is a rather unnatural construction, since it depends on a particular choice of basis. However, twirling 
with respect to local unitaries allows one to remove this artificiality and to generate a set of invariants. It turns out 
that these cumulant-based invariants are already known in the literature, though their relationship to cumulants is 
apparently not recognised. They account for five of the six independent 3-qubit invariants and eleven of the nineteen 
4-qubit invariants given in (l7| . Asymptotically, there are 0(2") of them whereas the total number of invariants in an 
n-party states is 0(2 n+1 ). There are also other families of invariants, including a family of 4th degree polynomials that 
accounts for about a quarter of all invariants. Breaking down the invariants into (overlapping) hierarchical families 
gives a perspective on their structure and, potentially, on the entanglement of states. 

There is alternative way of deriving invariants from cumulants (see Section IVlH A[) . due to Zhou et al [27| ■ Despite 
certain formal similarities, these seem not to have a simple functional relationship to our invariants. The entanglement 
space £ can be explored in other ways. Kraus [l3l [Llj showed how to reduce multipartite qubit states to a standard 
form that is invariant under local unitaries. An alternative, geometric procedure that applies to local spaces of 
arbitrary dimension (not just qubits) has also been proposed |22|. These methods enable one to determine if two 
states belong to the same orbit; invariants give a parametrisation of orbits and can be regarded as entanglement 
measures. All these methods give complementary insights into the structure of £ . 
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II. THE ALGEBRA OF MULTI-PARTITE STATES 

Let A^ be the the commutative algebra over C with generators e^, i = 1, . . . , n satisfying ef = 0. An element ip 
of An has the form 

where < ik < d — 1 and the a's are complex coefficients. Then ip = a + r, where a = ao...o, and r consists of at most 
d" — 1 terms and satisfies r n(d-i)+i = o. Any analytic function / : An ~> A^? can be expanded in a Taylor series 
in r: 

/W = /(a) + /'(a)r + /"(a)r 2 /2! + ... (1) 
This series is finite, because of the nilpotency of r, and is thus well defined. For instance, in An any ip can be written 

V' = a oo + aioei + a ie 2 + aneie2 (2) 

and if aoo 7^ 

log-0 = log(a 00 + r) = loga o + log(l + r/a 00 ) = loga 00 + r/a Q0 - r 2 /(2al ), 

= log a o + (aio/aoo)ei + (a i/a o)e 2 + («u/aoo - a w a i/al )eie 2 - (3) 

Similarly, there is a finite polynomial for the algebra inverse; e.g. 

V^ 1 = a oo _ ( a io/aoo) e i - (aoi/aoo) e 2 ~ (an/aoo _ 2ai aoi/ao )eie 2 ; (4) 
and for other functions, such as exp. These functions have all the expected properties, e.g. 

tin/,- 1 = 1, (5) 

log(^) =log(</>)+log(</>), (6) 

exp(V> + 4>) = exp(ip) exp(0), (7) 

exp(log)(V) = ip. (8) 

Now identify e^ 1 . . . e^™ with the n-qudit basis element \i\ . . . i n ). This sets up an isomorphism between elements ip 

of An and unnormalised n-qudit states \ip). For instance, ip = aoo + aioei + aoie 2 + aneie 2 in An can be identified 
with the two qubit state \ip) = aoo|00) + aio|10) + aoi|01) + an|ll). One can then carry across the structure of the 
algebra. For instance, if \<p) = &oo|00) + &io|10) + 601 101) + 6n|ll) we have the product 

\ip)\<j>) = a 00 6oo|00) + (aoo&io + ai 6 o)|10) + (a 00 &oi + a O i&oo)|01) 
+ (aoo»ii + a w b Q i + a i&io + an&oo)|ll)- 

The identity element in An , is |0...0), and the inverse, log and exponential are carried over in the obvious way 

n 

from the corresponding functions in A^ ■ 

Suppose the n subsystems are divided into two sets S and T. We write 

IV) = MS ® |V>T, (9) 

to indicate that \ip) is separable with respect to these subsets, the order in which S and T appear in the tensor product 
not being necessarily related to the order of their indices (e.g. we might write \ip) = \ip)iz <8> W)2)- Then l^) can also 
be written in terms of the algebra product as \ip) = \ip)s\ip)T, where \ip)s is identified with an element of the algebra 
that only uses those with i <E S, and \ip)r using only those e 2 ; with i G T. From (|6|) we have 

log \ip) =logh/>)s + log|#r. (10) 

thus linearising the tensor product when the log is defined, i.e. when the constant coefficients in the algebra do not 
vanish. 
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The coefficients c il .. An of log |0) = c ii-..i n \ii ■ ■ ■ in) are called cumulants. The cumulant corresponding to c^...^ is 
defined for classical random variables as the coefficient of A"..AJf in log(e^= A ^ x i) |2l| . This follows by identifying 
(X\ x . . . X n n ) with a^...^ and taking ao...o = 1- For instance, the equivalent of 

en = (anaoo - aioaoi)/aoo> ( n ) 

which is the coefficient of e±e 2 in (|3|), is the classical second degree cumulant (XiX 2 ) — (Xi)(X 2 ). The algebra Ah. 
can also be identified with the "moment algebra" in pj by associating to the term ce^ . . . ej n in .An the map that 
assigns to the integers {ii,...i n } the value c. 



III. MULTIPARTITE SEPARABILITY 



Write the state space for an n-party state as (C d )i <g> . . . <g> (C d )„. If 7r = {tti,tt 2 , ■ . . , 7Tfc} is a partition of n, we say 
)) is 7r-separable if we can write 



(12) 



i=i 



where each 10)^ is a state in the subspace ®j g7r . (C d )j. As we have seen, we can also write this using the algebra 
product as 

10} = 10)^...!^. (13) 

From ([6]) we get 

log |0) =log\(t>} ni +...+log\(f>) nk . (14) 

This immediately gives a characterisation of multipartite separability. Let us say that a set of indices i± . . .i n splits 
the partition tt if there are non-zero indices ij in more than one subset of tt. 

Theorem III.l. An n-party state \<p) with ao...o ^ is it-separable if and only if c^...^ = whenever ,i n splits 

TT. 

Proof. Necessity follows from the fact that the cumulants with indices splitting tt are absent from the expansion ((TJJ 
of log |0). Sufficiency follows by noting that, if these cumulants are zero, we can write log |0) in the form 

lo g|0)=X] X] Cilia. ..*JW2 ■ ■ -in) 

k {ij = l =^ jeir k } 

and exponentiating this shows |0) to be 7r-separable. □ 

The condition ao...o 7^ reflects the special role played by |0 . . .0) as the identity in the algebra. We shall shortly 
give a version of this theorem (|IV.2|) which does not have this unpleasant restriction. 

It should be emphasised that one can easily write down algebraic conditions for a pure state to be 7r-separable. 
However, the characterisation of Theorem IIII.ll will turn out to provide a useful starting point for making qubit 
invariants. It is also economical, in the sense that the vanishing of the c's gives the right number of equations to 
define the subspace of 7r-separable normalised states. Indeed, the (real parameter) dimension of this subspace is 

the expression in brackets counting the real and imaginary parts of each coefficient of {ij)}^, with 2 subtracted for 
normalisation and phase invariance. On the other hand, the number N v of index sets that split tt is 



this being the total number of non-empty subsets of {1, . . . ,n} minus those that do not split tt, i.e. those where the 
1-indices lie wholly within some 7Tj. But the total dimension of normalised states is 2d" — 2 and, as each equation 
Ci 1 ...i„ = contributes two constraints from the vanishing of its real and imaginary parts, we require 

dir = (2d n - 2) - 2N % , 



which is readily seen to hold. 
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IV. INVARIANTS FOR N-QUBIT STATES 

From now on we restrict attention to qubit states. We are interested in polynomial invariants of the action of local 
unitaries, which we take to be the group SU (2)™ x U(l). Thus, given a state l^) = i a^...^ \ii ■■ ■ i n ), we seek 
real-valued polynomials m the coefficients ai 1 ...i n and their complex conjugates di 1 ...i n satisfying 

I(g\i>)) = m)), (15) 

for any g £ SU(2) n x U(l). Equivalently, we can express the group action on the ith copy of SU(2) by 

Pi(9)aj 1 ...j i ...j n = Y^9j i k i aj 1 ...k i ...j n , (16) 

ft* 

and extend this to any monomial m — Hg=i a j q n by Pi(9) m = n {Pi(g) a j q i---j g „)■ We then require the poly- 
nomial / to be invariant under pi(g) for all i and g € SU(2) as well as invariant under phase changes introduced by 
U(l). 

The cumulant Ci 1 ...i n can be made into a polynomial di t ...i n by putting 

where 6 is the the number of l's in the set i\ . . .i n . This is a homogeneous polynomial of degree 6 in the a's. For 
instance, from pip we have d\\ = a^cn = anaoo — aioaoi- In general, we have the following formula for di x ...i n ■ Let 
S = {k\ik = 1}, so \S\ = is the degree of d^...^. If A C 5 let denote a^...^ where j^, = 1 if A; £ A and j'/. = 
otherwise. Then 

M 

di I ...*n = E(- 1 ) w " 1 (w- 1 ) !a !" |,r| II a '*' ( 18 ) 

7T 2 — 1 

where the sum is over all partitions tt = {tti, . . . , 7T|u-|} of 5, | vr | being the number of subsets in the partition. 
The action of local unitaries on c?n is given by 

Pi(g)d u = Ad lu (19) 

for i = 1, 2, where A = det g. To remove the dependence on the phase A, i.e. to get invariance under the action of 
U (1), we take 

/n = Mill 2 - (20) 

This gives us a local unitary invariant, and there is just one such invariant for normalised 2-qubit states. The 
transformation of cumulants for three or more qubits is more complicated (Theorem IIV.6[) . However, we can always 
get an invariant by integrating the squared modulus of di 1 ...i„ over the group SU(2) n , and this prompts the following: 

Definition IV. 1. 

Ii!...i n = 7n,e \pi(gi)---pn{g n )di 1 ...i n \ 2 dgi...dg n . (21) 

JSU(2)™ 

Here the integral is over the Haar measure, and the constant {actor j n< g = (8 + l) n ~ e (6 — l) 6 is introduced for later 
convenience. 

We now explore the properties of these invariants. First, we note that we can give a more satisfying, basis inde- 
pendent, version of Thcorcm lIII.il 

Theorem IV. 2 (Separability criterion). An n- qubit state \ip) is tt- separable if and only if 1^...^ = whenever i\ . . . i n 

Splits 7T. 

Proof. Suppose is 7r-scparable, so \tfj) = (Qlip}^- If a o...o 7^ 0, Theorem IIII.ll says that Ci 1 ...i n = and hence 
di 1 ...i n = 0. If a ...o = 0, define \ip) x = 0(1^}^ + ac|0 . . .0)^,). If x ^ 0, ao...o 0> so ^ti...»„ = f° r this state. 
By continuity, di t .,.i n = for x = 0, i.e. for the original \ip). Thus T% x ...% n = 0. Conversely, if = 0, then 

di 1 ..A n {g\4 ) )) = for all g, and for some g we must have ao...o 7^ 0, so Theorem IIII.ll can be applied. □ 
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Next wc define the action of elements g £ SU{2) on our cumulant polynomials. We will need some notation. 
Definition IV. 3. Define 5, by 

S.'tl ...:: .../ — aj 1 ...i i ...j n , (22) 
Sia h ... u ...i n =0. (23) 

Now, given a monomial m = Y\q=i a u q ...i n q °f degree 9, define TZ^f. to be the coefficient of x k in rj(i + x S^a^ q ...i n q ■ 
Extend this definition by linearity to any homogeneous polynomial of degree 9. 

For instance, for 9 = 3 we have 

TZ lA m = (5,11 + 15 t l + llSj) m. 

Note that, by virtue of the symmetry, this definition does not depend on the order of the a's in m. We can think of 
lZi t k as being like a raising operator (hence the 'R'). For example 

7^-3,0^110 = iiio a ooo — aioo°oiO: (24) 
72-3,1^110 = a m a ooo + ano°ooi — aioi°oio — aioo a oii: (25) 
72-3,2^110 = am a ooi — aioi a oii: (26) 

Lemma IV. 4. If li = 1 and at least one other index in di^...i n is 1, then TZi^ed^...^ = 1Zifl—idi- L ...i n = 0. 

Proof. Writing d for d; 1 ...z„, lZi : gd = because we cannot add 9 l's to the 9—1 a's that originally had zeros in 
position i. T^g-id has l's at position i in every a, so we can write lZi_g~id = d(|?/>)), where is a state with 

a>h...0i...l„ = fflJi...ij..J n - This means that factorises as ® \ip)j^-i. We now invoke Theorem IIV. 21 since 

■fji, . . . , l n } splits the partition {i}, {AT — i}. □ 

Lemma IV. 5. If we write Lid for the result of replacing all the 1 's in position i in d by 0, and if U = 1, then Lid = 0. 

Proof. The same argument as in the proof above shows that Lid = d(\ip)) where factorises. □ 

J . Assume 

that the index set l\, . . . , l n contains at least two 1 's. Then if U = 0, 

e 

Mg)d = Y,u e - k v k U hk d, (27) 

k=0 

and if li = 1, 

6-2 

Pl {g)d = A J2 u e - k - 2 v k K hk d, (28) 

fc=0 

where d stands for d^ ... j B . 

Proof. When U = and m = ITg^i a i 91 ...j»„ . 

= \\ («Oj e ,i...o < ...j 9 ,„ + uOj,,i...ii...i ll ») , (29) 

9=1 

and it is clear that the coefficient of u e is the original monomial m, and the coefficient of u e ~ k v k is TZi^d. 
When Zi = 1, 

Pi(ff)m = (wa jpA ... 0i ...j Pin + za 3pA ...i z ... Jvtn ) J| (™ J , il ...o,... J ,,„ + va jq l ... u ... jq n ) , (30) 

where a,j 1 j n is the unique a in m that has a 1 at position i. The coefficient of zv 8-1 is TZ^g—id, which is zero by 

Lemma TlV.41 and for k < 9 — 2 the coefficient of u 6 ~ k ~ 2 v k (uz) is cPi^d. The coefficient of u w is Ld, which is zero 
by Lemma riV.5[ and for k > 0, the coefficient of u e ~ k ~ 2 v k (vw) is TZi^+\{Ld) — TZi^d = — 72-j fed. Since A = uz — wit;, 
the theorem follows. □ 
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The polynomials TZi.kd appearing in equations (|27|) and (|28l) form an irreducible representation for SU{2). To see 
this, let h = ( ^ ^ ^ . Then, for = 0, applying to (|2"T)) gives 

^ u e " V ( Pi {h)1l itk d) = Pi(h) Pi (g)d = Pi {gh) = + vc) e - k (ub + vdfTl^d, (31) 

fe=0 fc=0 

and equating coefficients of u e ~ k v k in the left- and right-hand sides of (|3"Tj) gives the action of h at position i on all 
the polynomials TZi^d and thus defines the representation matrix, which can easily be recognised as the symmetric 
representation with Young diagram (9). When ij = 1, the same argument applied to (|28|) shows that we obtain the 
representation with Young diagram (9 — 1,1). Classically, cumulants were called "half-invariants" [f| Hf| because, if 
c is the classical joint cumulant in the random variables X\ . . . X n , then c is mapped to z n c when Xi is transformed 

by the affinc map Xi zXi + w. The equivalent of an affinc map in our setting is the group of matrices 



with x^O. From Theorem IIV.6I we see that the representation becomes 1-dimcnsional, sending d — > z 6 d. 
Theorem IV. 7 (Formula for invariants). // the index set i\, . . . , i n contains at least two 1 's, then 

2 

I, E ll ■■■:)''■ ' 

ki,...,k n p—1 



1 



Here, if i p = 0, k p ranges from to 9, and a° k = ( fc ) . If i p = 1, k p ranges from to 9 — 2, 



-2\ 



Proof. Combining (|21[) and Theorem IIV. 61 we get terms from f \pi(g)d\ 2 dg such as f \u\ 2i - e fc ' \v\ 2k dg\7Zi,kd\ 2 , and the 
integral can be calculated by using Schur's lemma: 

{ P+ P ^ I |u|2P|v|29d5 = { ^ J 9® {p+q) \0)(0\® ip+q H9T ip+q) dg\^ q ) (33) 

= dimSymP+^C 2 )- 1 ^,^^™!^.,) (34) 
= dimSym p+9 (C 2 r 1 (35) 
= (p + q+l)- 1 . (36) 

Here Ps y m is the projector onto the symmetric representation, and \ipp,q) is the normalised weight vector 
S CT IV(i) • • • V(p+g)); with ii = ... = i p =0, i p +i = ... = i p + q = 1, and with the sum taken over all 
permutations in S p + q . There are n — 9 indices i p that are zero, where (p + q + l) = 6 + l, and 9 indices i p that are 1, 
where {p + q + 1) = 9 — 1. These terms therefore cancel the constant 7^0 = (9 + l) n ~ e {9 — l) e in (|2"Tj) , and the a^'s 
come from the factor ( p+q ) in ([33|. 

There are also cross-terms in J \pi(g)d\ 2 dg, but these have the form f u^ e ~ k ^v k vS e ~^v^ dg (iZi^dlZi.jd) , and 

u («-*W«-jV dfl = (^_ fc , fc |^-«), (37) 



which vanishes for j ^ k since distinct weight vectors are orthogonal. □ 
As examples, consider the case of three qubits. From (|32|) we find 



hw = \K 3 ,od 110 \ 2 + i|^ 3 ,irf 110 | 2 + |^ 3 , 2 diio| 2 , (38) 

which is a polynomial invariant of degree 4. The terms are given explicitly by (|M|) . (|2"5"|) and (|2"B"1) . There are two 
other 4th degree invariants, iioi and Jon obtained by permuting the indices. We also have 

1 

hu = l^i^2j^ 3 ,fed m | 2 , (39) 

i,j,k=0 
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which is of degree 6. The theorem excludes the case where there is a single 1 in the index set. For instance, we have 
dioo = »ioo- If we apply Definition 1211 with 7 n l = 2", we find /ioo = (fp\fl>), and this result is independent of the 
ordering of the indices. This holds for any n, i.e. iio™- 1 = where fc stands for a string of k O's, and again the 

ordering of indices is immaterial. This invariant is often included in the list, assuming that states are not normalised. 

With this assumption, there are altogether six algebraically independent invariants for three qubits [HI, [23| . A 
further invariant is needed in order to have a Hilbert basis |9|, which has the property that every invariant can be 
written as a polynomial in the basis elements. However, the seventh invariant can be expressed in terms of Sudbery's 
six if one is allowed a square root, and from the point of view of a quantum information theorist rather than an 
algebraist, a complete, algebraically independent set serves to characterise the orbits. 

Following (23[, five of Sudbery's set of six are the second-degree polynomial J\ = (tpl'ip), the three fourth-degree 
polynomials J 2 = trp§, J 3 = trp^, J± = trpf and the sixth-degree polynomial J 5 = 3tr [(p 1 <g> p-i)p\i\ — tr(pf) — tr(p|)) 
where p\ = tr23 1 V) (V'l : P2 = t r i3 1 ^) | , pz = trialV'KV'li and p\2 = trz\tp){^j\- A straightforward calculation shows 
that our invariants can then be expressed in terms of these polynomials by 





— /010 


= /001 = Jl; 


4/no 


= J^ 


~ J2 — J3 ~ J±i 


4/ioi 


= jN 


' J3 — J2 — Ji, 


4/on 




- J4 — J2 — J3, 


6/m 


= 5Jf 


- 3Ji ( J 2 + J3 



Ji) +4J 5 . 

Alternatively, they can be identified with five of the list given in fl7| : /001 is their ^4m, /no is B002, and /m is Cm. 

Next we consider the question of independence of the invariants. We can quickly eliminate certain types of functional 
dependence by using the separability criterion, Theorem IIV.2I For instance, we cannot have /no = /(/in, /101, /on) 
for any function / because all the invariants on the right-hand side are zero for any {12}{3}-separable state, whereas 
/no can take a range of values according to the choice of the state. However, we cannot use the separability criterion 
to rule out a relation of the form /m = /(/on, /ioi ; /no)- say. Our next aim is to prove that the invariants are in 
fact algebraically independent, i.e. that there can be no non-trivial polynomial relation between them. 



^^ il<1 ... ik ji^,...ii", =0, (40) 



Theorem IV. 8 (Algebraic independence). A polynomial relationship between the invariants 1^ 

En- T jl P k 

U 8l,l—«fc,„- t 8l>1 ... lli „ ■ ' ■ - t ifc,l...'ifc,„ 

can only hold if each 1 ...i k „ =0. 

Proof. Consider the neighbourhood of the fully separable state |0. . .0), where the ^...j^'s with not all ij = are 
small. From (|3"2")l and (fT5)) we see that the expansion of the invariants in lowest degree terms in these small variables 
only has contributions from the term aji...j„ao~o m ^ji - in; and from the corresponding term where operators TZ^i 
applied to positions where the index ji = 0. Thus, assuming that the indices ii,...,i n include at least two l's, we 
have 

1 1 

u^uE- E i" >: ' ) •'•< •••< • , • ( 41 ) 

n 1 =0 »t„_ 9 =0 

where Xj 1 ...j n = \dj 1 ...j n | 2 , and t\ . . .t n -g denote the positions of indices in i\, . . . ,i n that are zero. We also have the 
special case 

1 1 
/10...0 = E • E ■■■*»• 

i 1= i„=0 

Let us write I\ = /10...0 — ^o.-.O! and x\ = X10...0 + 2^010. ..0 + • ■ ■ a;o...oi- Then the lowest degree expansion defines an 
invertible map A. For instance, for n = 3 this map is given by 





h 


/no 


/101 


/on 


/111 


Xi 


/I 














XlW 


1 


£000 











X101 


1 





£000 








Xoil 


1 








£000 





Xlll 


V 1 


2" x 000 


5^000 


J^OOO 


x ooo 
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This extends to a map A poly from polynomials in the Ps to polynomials in the x's. Since A is invcrtible, the same is 
true of A poly . If we now take the lowest degree terms in (|40p . they map to a polynomial in the cc's whose coefficients 
must be zero, since the x's are independent variables. Applying the inverse of A poly , we deduce that the corresponding 
coefficients in P0|) are zero. Looking at the next highest degree terms in (fH))) , we again deduce that their coefficients 
are zero. And so on for the whole /-polynomial. □ 

We can therefore conclude that we get 2™ — n independent invariants from cumulants, namely one invariant, the 
squared-amplitude, from the first-degree cumulants, and 2™ — n — 1 from the second and higher degree cumulants. On 
the other hand, the total number of independent invariants [|| for n > 3 is iVj = 2" +1 — (3n+ 1). Thus asymptotically 
there are ATj/2 cumulant-bascd invariants. 



V. THE HIERARCHICAL STRUCTURE OF INVARIANTS 

The 3-qubit invariants ioiii-^ioi and I no are closely related to the 2-qubit invariant In- We obtain Ino from In 
by adding a '0' index in the third position to In and twirling. More generally, we wish to make an (n + l)-qubit 
invariant from an n-qubit invariant J. Phase invariancc requires every term in J to be a product of an equal number 
of a's and a's. Given a monomial m = Y[ a ji---jn®k 1 ...k„ in J, of degree 9 in the a's and a's, we define 



m = (9 + 



!) / /',':.'/;] [".. ...o ...j a kl ...o .../. dg. (42) 

JSU(2) 



We refer to this process as lifting and use a subscript '0' to denote a lifted index. Table U shows lifts of In and 7m 
up to five qubits. We say that an invariant together with its lifts constitute a family. 

Lifting can also be understood in terms of tracing-out operations. Let us write a mixed state of n qubits as 

n 

p= E <::t(g)>><>: ( 43 ) 

il,...»niJl,— On r=l 

The following rules, applied to each monomial, enable one to interconvert between a pure state invariant J of degree 
9 in the a's and a's, and a mixed state invariant J of degree 9 in the coefficients of (|13"|) : 



1 c 

•/ II «*.!....« II ^..-iS ^E II 4,^" J " " ( 44 ) 

r—1 s—1 <7$LSo r— 1 

j -> j • li — ► II ^.-.iK II %T,-,is (45) 



r— 1 r—1 r—1 

Proposition V.l. If pi(g)J = J then pi(g)J Pi{gy = J- 
Proof. The following diagram commutes 

J J 

I I 

J J 

The following is immediate: 
Proposition V.2. 

Jdvxvi) = Jim- 

Proposition V.3. If is an n-qubit state and J Pl ... Pk a polynomial invariant, then 

J pi ...p k 0"-«(\^)) = Jpi... Pk (^n-k\lp)W), 

where n ~ k means a string of (n — k) '0' indices, and tr„_fc means tracing out the last (n — k) systems. 



□ 
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Proof. Let us sec how this works in a simple case. The generalisation is then straightforward. So consider the 
monomial m = ai li2 aj 1 j 2 ak 1 k 2 di 1 i 2 01 a 4th degree 2-qubit invariant J piP2 . Under the map (|44[) m becomes m = 

i (af\ k2 a 1 ^ +al ll J i a^A. The monomial m lifts, by gl, to 

mo = (0 + 1) y P3(g)ai 1 i 2 oa jl j 2 oa kl k 2 oai 1 i 2 odg (46) 

in ■/ P ip 2 o(|'0)) f° r a 3-qubit state |V>), and (0 + 1) = 3. We therefore wish to show that too = to, where the coefficients 
a iii 2 2 ' e * c -' m ^ come from tr3|^)(^|. This gives 

m = 2 ( a u*2oafci/c 2 o + ai 1 i 2 ia,k 1 k 2 i) ( a jij20^hho + a jij 2 i^hl 2 i) (47) 

+ 2 ( a ii»2oa/iZ 2 o + a ili2l ai 1 i 2 i) (a,j 1 j 2 oak 1 k 2 o + a jij 3 iflfeife 2 i) ; 

and from (l46l) we have 



to = J (ua lll2 o + va ili2 i) (uaj lj2 o + va Jlj2l ) (uak 1 k 2 o + vak 1 k 2 i) {uai^o + vai^i) .dg (48) 

Comparing (|4"T)) and (|4"3j) . and using ([33]), ([35|) and (|37|) . we see that too = to. This argument is unchanged when J 
is an n-qubit invariants (we just permute notationally more cumbersome blocks of indices). When the degree, 28, is 
arbitrary, a term in to that has p O's in the a's in the lifted index position occurs p\{9 —p)\ times in the generalisation 
of (|47|) . With the normlising factor 1/0! from (|44]) . we obtain the coefficient J \u\ 2p \v\ 2e ~ 2p of the corresponding term 
in mo, as given by (|35]). □ 

Proposition IV.31 provides an alternative definition of lifting via tracing-out of subsystems. Combining this Propo- 
sition with Proposition IV. 21 we get 

Corollary V.4. If\ip) = \fj,) k (8 \f) n -k then 

Jh...^-"^)) = Jh-iAlfj))- 

The third, equivalent, definition of lifting comes from a technique called transvection, invented by Cayley in the 
19th century heyday of invariant theory. Transvection is a useful device for generating invariants; understanding it 
will enable us to interpret the 4-qubit invariants given in [TtJ in the language used here. 

The fundamental form for a n-qubit state is the polynomial / in the a's and the variables x^ , , for 1 < j < n 
given by 

/= £ a,..,.^...^ (49) 

If we let g € SU(n) act on the ith index of a's by the usual transpose action (|16[) and upon the #W >g via the inverse, 
g', then one easily checks that p(g)f = f ■ More generally, a covariant of weight q is a polynomial p in the a's and x's 
satisfying 

p i (g)p = A q i p. (50) 



Given two covariants, p and g, we define (p, g), by (/i(a)|^(a)) = ji(a)v(a) for expressions in the a's, and, for each i, 

{(x^rix^^Yix^Y) = W, (51) 

This is sometimes called the derivative inner product because we obtain it by setting x^ on the lefthand side to 

d/dx^ and applying these derivatives to (unchanged) x's on the righthand side. 

If p is a covariant, then (p,p) is an invariant, so any means of generating covariants also supplies us with invariants. 
Transvection is just such a means. Given two covariants, p(Xj ), q{y^), define the transvectant by 

(p, q y^ = <>, ...!>..:>,; | where 0* = - 44- (52) 

y~* x ox oy\ ax\ oy 

The vertical bar indicates that, after applying the differential operators Qj, we change the y's to x's, so (p, q) %1 " Ak 
is a polynomial in a's and x's. A classical theorem fl9j asserts that, for any binary indices i\ ■ ■ - i k , (p,q) ll "' lk is a 
covariant if p and q are. Starting with the fundamental form, we can build up a wealth of covariants p and derive 
invariants (p\p) from them (see Table |H]) . 
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TABLE I: Some 2, 3, 4, and 5-qubit invariants related by lifting. 



2 qubits 


3 qubits 


4 qubits 


5 qubits 


ho 


hoo 


hooo 


/lUOOO 


hi 


hw, hoi, Ion 


hum, how, /0110, hool, ^oioi, /0011 


^11000,^10100, etc. (10) 






Gim 


Gi 1 no, G11 101 , Gi 1011 , Giom , Gomi 




/in 


hno, hioi, hou, Zdhi 


/moo, /11010, etc. (10) 




H222 


H.2220 , H2202 , #2022 , #0222 


#22200, #22020, etc. (10) 



TABLE II: Sixteen of the nineteen four-qubit invariants in transvectant notation. 



Invariants 


Corresponding covariants 


number 


degree 


hooo 


/ 


1 


2 


Gim 


(/,/) mi 


1 


4 


J1100 


(/,/) 1100 


6 


4 


hno 


^1100)0010 


4 


6 


#2220 


(/,(/,(/,/) 110 T 010 ) 1110 


4 
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Example V.5. For two-qubit states, f = Y^ a ij x i^ x< f^ ■ Take p = (/, /) u . Then p = d\i and (p\p) = \dn\ 2 = In. 
For three-qubit states, f = Yl a ijk x i x j x k ■ Take tno = (/, /) 110 - Applying \5S\) we get 

Hio = (/,/) n ° = rfno(4 3) ) 2 +^3,idno(4 3) 4 3) ) +^3, 2 diio(a;i 3) ) 2 . 
Using the derivative formula for the inner product \51\) we find that (tiio|<-iio) = 4/no- 

More generally, we have the following result: 
Theorem V.6. Let 

tl , „- fc = (/,(/,... (/,/) iio - o r-°...) o - oio "- fc . 

Then 

where £ = 4((fc - 2)!) fc (fc!)"- fc . 

Proof. Consider first the terms in u-^Qn-k where the subscript in every x is O.The first transvectant step, (/, /) 110 '°, 
yields terms 

(aiii 3 ...i n a o h ... jn - a 10 i 3 .. An a 01h ... jn ) xfjxfj . ..x^xfj. (53) 
If k = 2, the restriction to ccn/s means that we get 

I J ( (3) (")A 2 

t ll(0"- 2 )ko — «11(0™- 2 ) 1^0 ---Xq j . 

If k > 2, at the next transvectant step we set the x's in (|53[) to y's, multiply by the fundamental form / and apply 
f^3 to get 



E 



(i) 

a fel ...fc„4/ 



(n) 



E 



(axii 3 ...t„oooj 3 . 



N (3) (3) 

aioi 3 ...i„aoy 3 ...in)%3 % 



(n) (n) 
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If we are restricted to xq j s, we must have k\ = ki — since no further £1 operations arc applied in these index positions 
and so these x's will be unchanged. Only certain sets of indices are consistent with a y remaining after applying tt 3 

to x^y^yfj; namely (1) k 3 = 0, i 3 = 1, ja = 0; (2) fc 3 = 0, i 3 = 0, j 3 = 1; (3) k 3 = 1, * a = 0, j 3 - 0. The result of 
this operation is of the form 

X^(rv- • , -I- R- ■ , W (1 M 2 M 3 M 4) t. (4 M 4) >'t ( "'tW 

where otijki @i,j,k are terms in the a's with compound indices i = {i^ . . . i„}, etc., and onjk comes from the conditions 
(1) and (2) above on index sets: 

<*i,j,k = aooofc 4 ...fc„7?-3,i [aiioi 4 ...i„ a oooj 4 ...j„ — aoioi 4 ...i ri aiooj i ...j n ] , 
whereas from condition (3) we get 

Pi,j,k = aooifei...fc 4 [ a iioi 4 ...i„aoooj 4 ...j„ — aoioj 4 ...i„aiooj 4 ...j„] ■ 

If k = 3 this simplifies to 

iiii(0"-3)Uo = Koo(o™-3)^3,idiio(0"-3) - 2a oi(o»-3)diio(o»-3)] (4^ 'x { q ] '4, 3) ) (x { i} . . .4"') ■ (54) 

Using (|18p . a straightforward calculation shows that 

diii(o»- 3 ) = aooo(o™- 3 )^3,i^iio(o™- 3 ) — 2a 00 i(o™- 3 )rfiio(0' 1 - 3 )- (55) 
Repeating the above argument, we have 

Llh n-k\ Xo = dlk n-k (xq ^ . . . Xq M (xq + ' . . . Xq M , (56) 

and the generalisation of (|55j) is 

dlkQn-k = aonlZk,ldlk-lQn-k + l — Idftk - 1 J^Qii - k + 1 rfj fe - 1 Qn .— k + 1 . (57) 

This last equation has a straightforward interpretation. When evaluating f(ip) by the Taylor series ([1}, the coefficient 
of, say, ei can be obtained by differentiating ^-/(?/0 and setting = 0, for all i. Writing ip = a + r, where a = a , 
we find 

— /(a + r)| ei=0 = / (a + r) — r| £i=0 = / (a)aio...o = / (a)fti , x a. 
aei aei 

We can interpret the last expression above as the formal derivative of using the raising operator 7^1,1, and 

similarly the coefficient of any product . . . ei q is the result of formal derivatives by TZi 1 ,i ■ ■ ■ TZ-i q .i- This can indeed 
be taken as the definition of the expansion of f(ip), as in For the log function, the coefficient of ei . . . e^-i is 
C^fc — 1 Qn— fe+1 , and the coefficient of e\ . . . e^, namely c^Qn-k is obtained by applying IZk.i to c Differentiating 
log(V0 an d using c lqQn - q = d 1 q n- q (aon)~ q gives ([57]). 

From (|56"|) and the definition of the inner product (|5ip we find that 

((ti*o»»-*Uo) I ( t i fc o™- fc Uo)) = €Mi fc 0"-- fc | 2 j 

where £ is the constant given in the Proposition. With the restriction to xo's we therefore get, up to the factor £, 
the term in the formula for Iik n-k (Theorem IIV. 71) where k p = for all p. To complete the proof, one observes that, 

(i) 

allowing k x\ s introduces k l's into the a's at position i, and is equivalent to applying lZi pi k p to d^Qn-k. The values 
of the coefficients a£ are given by the derivative inner product. 

□ 

This enables us to recognise some of the four-qubit invariants in [l7| . Up to a constant factor, our iioocb ^noo 
and /mo correspond to their Ami, (-B0022I-B0022) and (C1113IC1113), respectively. We use different letters for the 
invariants, and our subscripts indicate the total number of l's at a given position in successive transvection operations; 
see Table HU 

We now come to the third way of defining the lift. Suppose that a covariant pi 1 ...i n is derived by some sequence of 
transvectant operations. Define its ith lift p; 1 ...o i ..j n by adding an index position in the ith position in the ground 
form, and applying the same transvectant operations, but with an '0' added to the transvectant indices in the ith 
position. 
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Proposition V.7. If Pi t ...i n = (pi 1 ...i n \pi 1 ...i n ) is the invariant derived from the covariant pi lm ..i„, then the ith lift of 
Ph. ..l n ^ given by /',..,, , = (pi 1 ... 0i ...i n \pi 1 ...o l ...i n ) ■ 

Proof. Because there is a at position i in the transvectant indices, fij is never applied during the transvection 
operations. This means that we get all possible products of x ^ and xf 1 , and the terms with k 's correspond to 
products of a's with k l's in index position 1. □ 

As an example, consider the invariant of highest degree for 3-qubits in [l7j . It is given (in our notation) by 
-H222 = (/J222 ^222), where 

^222 = (/,(/,(/,/) 11 °) 001 ) 111 . (58) 

Equivalently, #222 = |D et(h />))| 2 , where DetQip)) = aijkaiiji m anpk>a n ipi m itii'£jj>tkk>tmm>£nn>£p P > is the hyperdetermi- 
nant [l8|]. JJ222 is Iq m [IH, an d i s closely related to the 3-tangle r = 2|Det(|V>})|- From the lift of /1222 at 
position 4 is /i 2220 = (/,(/,(/,/) 1100 ) 0010 )" 10 , and therefore by Proposition E3 #2220 = (kraolkroo) is the lift at 
position 4 of H 2 22- In the terminology of JT3, #222 is £>ooo and J/2220 is (JJ0004I J^oocw)- 

Another example is obtained by putting 31111 = (/, /) and setting Gnu = (<?mi|<7mi)- This 4th degree 
invariant can be written 

Gnu = aooooaim - (aioooaoixi + permutations ) + (onooaooii + permutations ) . 

In general, for each k, we add a new 2fc-party invariant G^k = (g^klg^k) , where g^ik = (/, f) 1 . Together with all 
its lifts, the G family comprises (2) + (") + (") + . . . = 2"~ 1 — 1 independent invariants of degree 4 for an n-qubit 



system. This coincides with the family Bd in [17 



VI. CONCLUSIONS 

We have arrived at a view of local unitary invariants of qubit states which can be summarised as follows: 
• Half of them (in the limit of large n) can be taken to be twirled cumulants. Thus for n = 3 the cumulants 



J100, J110 and /in correspond to Am, -B002 and Cm, respectively, in [17} , which, together with permutations, 
account for 5 of the 6 invariants; and for n = 4, /1000, Iiwo and Imo correspond to their Aim, (-B0022IS0022) 
and (C1113IC1113), respectively, which, with permutations, gives 11 of the total of 19. 



All the invariants come in families, related by a tracing operation (Proposition IV. 31) that I call a lift and 
indicate by a zero in the index string. The cumulant family is obtained this way, from I\, In, /m, etc., by 



lifting. A further example is the hypcrdctcrminant family: for n = 3, Dqoq in the notation in 17J is the 3- 
tangle (the modulus-squared of the hypcrdctcrminant), which I denote by H222, and their (-D0004I-D0004) and its 
permutations are lifted 3-tangles, #2220 e tc. in my notation. These families may overlap. For instance, Jim 
is not included in the list of primary invariants in [l7| . and is therefore algebraically dependent on the lifted 
3-tangles and others in their list. 

Many of the invariants are closely related to separability of states. Thus, the vanishing of members of the 
cumulant family can be used to characterise multipartite separability (see Theorem IIV.2|) . The hypcrdctcrmi- 
nant, in its guise as the 3-tangle, is of course also an entanglement measure. Unlike the cumulant invariants, 
the 3-tangle measures entanglement of mixed states, whereas the mixed state equivalent of cumulants, J, only 
measure correlation in mixed states. 

We can add a number of conjectures to this list: 

• Jn attains its maximum of 1/4 on Bell states, and Jno is maximised by <S) |0), with *5> a Bell state. We can 
regard this as an example of monogamy of entanglement 0, E^j with Jno detecting entanglement between the 
first two systems, which achieves its maximum when they are unentanglcd with the third system. I conjecture 
that JikQ is maximised by states of the form \fj,) <£) |0), where the fc-qubit state \fj,) maximises J x a . Equivalently, 
Jit (p) is maximised when p is pure. 

• For 3 qubits, the highest degree invariant in [Tt], HI] is #222, (equal to |r 2 , where r is the 3-tangle HI) is 
derived from the covariant im underlying Jm by applying what one could call the total transvectant operation 
^222 = (/j tui) , with an il operation at each index position. Can a highest degree invariant always be derived 
this way, by ft.2« = (/, ii") 1 "? In [l7|, the highest degree invariant for 4-qubits is defined by (/, (/, iom) 1011 ) 110< \ 
but we can replace this by iJ2222 to obtain an algebraically independent set (though not necessarily a Hilbcrt 
basis). I therefore conjecture that J^" constitute another family. 
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How far can the ideas here can be generalised beyond pure qubit states? The cumulant-based invariants can be 
applied to mixed states via the map (|44[) . However, they then constitute only an exponentially small fraction of the 
estimated 2 2n — 3n + 1 invariants for mixed qubit states [24| , and furthermore Theorem IIII.ll tells us only about 
correlation rather than mixed-state separability. The results also fail to generalise for pure states where the local 
dimension exceeds two. We can construct invariants, and Theorem IIII.ll holds, but the invariants are not algebraically 
independent. This is seen even for two qutrits, where we have four members of the cumulant family, namely In, I12, 
I21 and ^22, whereas there are only two independent invariants [Toj . Since four polynomial equations is the correct 
number to characterise separability, the simple relationship between invariants and separability cannot hold for d > 2. 
Nonetheless the basic concept of lifts and families still applies in all these wider contexts. 
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VIII. APPENDIX 

A. An alternative cumulant-based invariant 

There is a very different way of relating cumulants and invariants, due to Zhou et al. [27j . Given an n-party mixed 
state p, one defines its cumulant by analogy with (fL8|) as 



M 

(-1)1-1-^1 -!)!(£) P7ri , (59) 
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where p ni is the result of tracing out from p all systems apart from those with labels in 7Tj. For instance, for three 
systems 



The cumulant operator given by is not in general a state, but Zhou et al. propose M(p) = ^tr\p c \ as a measure of 
correlation of the mixed state p. It is manifestly invariant under local unitaries, and, because of the general property 
cumulants have of vanishing on products, M(p) = whenever p = ps ® pr- For pure states, this means it vanishes 
when states are separable. 

It therefore seems to have formal similarities to our cumulant-based invariants, and one can carry this further by 
defining, in line with Proposition IV. 3[ the lift of M to be M (tr\ip) We can in fact adopt parallel notation to the 
7's, writing, for a 3-qubit state for example, Mm(\ip)) = |tr| (\4>}(ip\) c |, M lw (\ip)) = itr| (tr 3 |i/ , )(V'l) c I, an d so on. 
Then Mi 1 ...i n (\tp)) = for any 7r-separable \tp) where {i\ . . . i n } splits n. Furthermore, for 3-qubit states Mm(\tp)) = 
is sufficient for separability of \4>) ([13, Theorem 3), and the same is true if Im(\ip)) = 0. 

These similarities prompt the question of whether there is a functional connection. Can one write M^...^ (|^>)) = 
F(Iii-..i n (\ip))), for some function F? For 2-qubit states, M\\ = In + y/Iu. However, there is only one 2-qubit 
invariant for normalised states, so a functional relationship here is unsurprising. For 3-qubit states of the form 
|V') = a|000) + fejlll) one finds 



Pc = P~{pl® P23 + P2 <E> Pl3 + P3 <8> Pl2) + 2/9i ® P2® p3- 




(60) 



whereas, for states of the form \<j>) 



a|100) +6|010) + c|001) 



(M U1 - '-ff - ij m 



= 0. 



(61) 



Since (pO)) and (|6 1 [) do not define the same function of im, Mm must depend on other invariants besides Jm. 



